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The CKM Matrix from Lattice QCD 

Paul B. Mackenzie 
Fermilab, Batavia, IL 60510, USA 

Lattice QCD plays an essential role in testing and determining the parameters of the 
CKM theory of flavor mixing and CP violation. Very high precisions are required for lat- 
tice calculations analysing CKM data; I discuss the prospects for achieving them. Lattice 
calculations will also play a role in investigating flavor mixing and CP violation beyond the 
Standard Model. 

§1. Introduction 

ON 

The Cabibbo-Kobayashi-Maskawa (CKM) matrix 1 )' 2 ) parametrizes the couplings 
between flavors of quarks under the weak interactions. Quarks are believed to be 
permanently confined within hadrons, and the dynamics of quarks and gluons in 
hadrons is nonperturbative. Lattice quantum chromo dynamics (QCD) is the only 
general method for performing nonperturbative calculations in QCD, so it is not sur- 
prising that many analyses of the physics of the CKM matrix require lattice QCD 
calculations. 

If the Standard Model of particle physics is solely responsible for flavor mixing, 
the CKM matrix will be unitary, and parametrizable with four parameters, such as 
those of the Wolfenstein parametrization, 3 ) A, A, p, and fj: 

1-^A 2 A A\ 3 {p-irj) 

r= I V cA V cs V rh 1^1 -A 1-±A 2 ,4A 2 I (!■!) 
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The rows and columns of unitary matrices are orthogonal, satisfying 

k k 
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The off-diagonal combinations vanish, and can be represented by triangles in the 
t\ complex plane. The best-known of these is 

> 

v ud v: b + v cd v; b + v td v t * b = 0. (i-3) 



If this equation is divided by V c dV*b an d parametrized as in the right side of Eq. ( 1-1 ), 
one obtains the familiar p-rj unitarity triangle of Fig. [T] If the CKM matrix is the full 
story of flavor physics and the matrix is unitary, measurements of p and rj should 
all be consistent, whether they are obtained via V ub , Vtd, or through some other 
constraints. Conflicting determinations of p and rj would constitute evidence for 
physics beyond the Standard Model. 

Lattice analyses of meson decays and mixings provide some of the most impor- 
tant ways of determining CKM matrix elements. Fig. [T] shows some of the most 
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Fig. 1. Bounds in the p-fj plane from a global fit. (CKMfftter. 4 ') Lattice QCD is responsible for 
the bounds from KK mixing (light green band) and from BB and B S B S mixing (yellow and 
orange bands). It helps with the bound from \ V u b\ (dark green band). 

important constraints on ~p and f]. In KK mixing (light green band), and in BB 
and B S B S mixing (yellow and orange bands), lattice QCD is the only first-principles 
way of determining the nonperturbative parameters that relate the quark scatter- 
ing amplitudes of the Standard Model Lagrangian to the meson mixing parameters 
observed by experiment. In meson decays, determinations of CKM elements from 
nonperturbative lattice calculations of leptonic and exclusive semileptonic decay am- 
plitudes provide an important complement to determination from inclusive decays 
via perturbative QCD. All of the CKM matrix elements except Vtb can be determined 
from one of these exclusive processes with lattice QCD. Table [TJshows the CKM ma- 
trix elements and the process to which each is most sensitive: meson leptonic decay 
constants, Jm, exclusive semileptonic decays, M\ — > M2I1S, and meson-antimeson 
mixing amplitudes, < M\M >. 

Lattice Quantum Chromodynamics is a tool for understanding strongly inter- 
acting QCD at low energies, where QCD perturbation theory fails. It is an essential 
tool for studying weak interactions in low-energy hadronic physics because the effects 
of weak interactions involving quarks are shrouded at low energies by the effects of 
the strong interactions. These strong-interaction effects must be understood quan- 
titatively before the weak interaction properties can be inferred. For example, when 
a meson decays weakly into leptons, the decay amplitude is proportional to both a 
CKM matrix element and a hadronic decay constant, which parametrizes the ampli- 
tude for the two valence quarks in a meson to interact at a point. 

In lattice QCD, quantum fields for quarks and gluons are defined on the sites and 
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links of a four-dimensional space-time lattice. The quantum fluctuations of the fields 
described by the field theory path integral are calculated with Monte Carlo methods. 
The physical theory is defined as the zero-lattice-spacing limit of the lattice theory. 
Because the number of degrees of freedom in the path integral becomes infinite in 
this limit, this limit is computationally expensive. 

In the last decade, due to improvements in methods, in algorithms, and in 
computers, lattice QCD calculations have become able to produce serious, first- 
principles results for many simple but important quantities. Prime among these for 
CKM physics are the decay constants, exclusive semileptonic decays, and meson- 
antimeson mixing amplitudes of stable mesons. Stable pseudoscalar mesons are 
among the most tractable quantities for current lattice methods. They have the 
smallest and best controlled uncertainties for statistical errors, finite volume errors, 
and other quantities. Since they also provide some of the most accurate CKM- 
related experimental results, they provide the most accurate determinations of the 
CKM matrix from lattice QCD. 

Current lattice calculations are done on computers that are a factor of 10 8 more 
powerful than the VAX ll/780s on which the first numerical lattice QCD calcula- 
tions were done. Remarkably, just as large a factor of improvement has come from 
improvements in algorithms, and an even larger factor has come from improved meth- 
ods. This has made it possible to abandon the tactic in early lattice calculations of 
ignoring the effects of quark-antiquark pairs (the "quenched" approximation). The 
most serious lattice calculations are now all unquenched. Lattice calculations con- 
tain errors due to discretization, to extrapolation to the chiral limit, to operator 
normalizations, etc. Estimating the expected precision in light of these uncertainties 
is a key element in serious lattice calculations. The errors due to discretization have 
been greatly ameliorated by the use of improved actions. There are a half dozen 
families of actions for discretizing fermion fields with widely disparate virtues and 
drawbacks. Which, if any, is best for any given purpose is still contentious. The 
good news for observers and customers of lattice calculations is that many impor- 
tant quantities are becoming available in several types of quark methods, so that 
outsiders can see for themselves how well they agree. 

Most quantities that have been calculated so far with lattice QCD agree well 
with nonlattice results to the expected precision. A few, however, do not, and 
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Table I. All of the CKM matrix elements except Vtb can be determined from a tractable lattice 
QCD calculation. 
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these are some of the more interesting cases in current lattice phenomenology. For 
example, lattice calculations of fo s , the leptonic decay constant of the D s meson, lie 
significantly below the experimental results, while lattice calculations of other decay 
constants are spot on. Determinations of CKM matrix elements from inclusive B 
decays differ by around two sigma from those with lattice calculations of the exclusive 
semileptonic decays B — ► D*lv and B — ► nlv. I will describe these results as well as 
the much larger body of results for which everything is very consistent. 

§2. Lattice Determinations of the CKM Matrix Elements 

Flavor physics was reviewed at Lattice 2008 by Lellouch 5 ) (for kaons) and by 
Gamiz 6 ) (for heavy flavor). I will often refer to their talks for the state of lattice 
CKM data, and will often use their averages rather than provide a complete review 
here. I will discuss only unquenched lattice results, mostly with three light flavors, 
sometimes with only two. I will discuss all of the CKM matrix elements, but I will 
go into more detail where there is something interesting to say, and discuss others 
more briefly. 

In this article, I am grouping together CKM matrix elements which can be 
derived from related physical quantities and calculated with related methods. For 
each set of elements, I will then discuss methods to determine them. In the next 
four subsections, I discuss: 

• V c d and V cs . These may be obtained from lattice calculations of leptonic and 
exclusive semileptonic decays of D and D s mesons. 

• V u b and V c b- These may be obtained from leptonic and exclusive semileptonic 
decays of B mesons. 

• V u d and V us . These are obtained from pion and kaon leptonic and semileptonic 
decays. 

• Vtd and Vts- These are obtained from the meson-antimeson mixings < B\B > 
and < B S \B S >. I also discuss < K\K > mixing in this section. 

These quantities constitute the core of lattice CKM phenomenology. 

2.1. V cd and V cs 

\V c d\ and \ V CS \ can be derived from leptonic and exclusive semileptonic decays of 
D and D s mesons. In both cases, the processes are related to each other by SU(3) 
flavor symmetry. 

2-1.1. f D ,f Ds 

fo s presents one of the very few disagreements between lattice calculations and 
other results. I will therefore discuss these calculations in detail. The HPQCD col- 
laboration has calculated the four decay constants /„-, fx, Id, and /a/' with an 
improved form of staggered fermions called "HISQ" quarks ( "highly improved stag- 
gered quarks" ) . 8 ) Staggered fermions have multiple poles in quark propagators which 
cause a doubling of quark flavors, and whose effects must be removed from physical 
calculations. Transitions between these multiple poles (or "tastes" in lattice jar- 
gon) cause significant discretization errors. The usual improved staggered fermions 
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(so-called "asqtad" fermions) remove these taste breaking effects at the one-gluon 
level. HISQ quarks remove two-gluon scattering effects as well. In addition, to apply 
HISQ to charm quarks, HPQCD has calculated the 0{ma p ) errors to a high order 
and removed them, ((ma) 4 turned out to suffice for sub-per cent precision.) These 
improvements allowed the calculation of the four decay constants to higher precision 
than ever before. They obtained : 7 ) 

U = 157(2)MeV (2-1) 
U/f K = 1.189(7) (2-2) 
f D = 207(4)MeV (2-3) 
f Ds = 241(3)MeV. (2-4) 

The first three agree with experiment to within the stated precisions, a few per cent, 
but fjj e presents a puzzle, as I now discuss. 

fir, Ik, Jd, and fo s have also been calculated by the European Twisted Mass 
Collaboration using twisted-mass fermions. Fermilab/MILC have calculated fo and 
fu a using asqtad light quarks. These groups obtain results compatible with HPQCD, 
but with less precision. ETM obtain 9 ) f Ds = 244(8) MeV and f D = 197(9) MeV. (I 
will return to the subject of pion and kaon decay constants in a later section.) Fermi- 
lab/MILC, using clover /Fermilab charm quarks and asqtad light quarks obtain 10 )' 11 ) 
f D = 207(11) MeV and f Ds = 249(11) MeV. These results are shown in Fig. [2} Fig.|] 
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Fig. 2. Theory and experiment for fn and }d s - Theory results from HPQCD, ' Fermi- 
lab/MILC, 10 '' n) and ETM. 9) Experimental results from CLEO 12 ' (f D ) and HFAG 13 ' (f Ds ). 

also shows the value of fn from CLEO-c, 12 ) and a recent preliminary HFAG world 
average, 13 ) fjj a = 263.9(6.7) MeV. The results for theory and experiment are nicely 
compatible for fjj (as they are far the vast majority of lattice calculations), but 
disagree significantly for fjj a . This discrepancy is the largest discrepancy that has 
arisen in lattice phenomenology. It may be resolved by the lattice or experiment 
moving or changing their uncertainties. The experimental average has come down 
in the last year, and the discrepancy was 3.8 sigma a year ago. A recent CLEO-c 
number is lower still, quoting 14 ) fu a = 259.5(7.3) MeV, so it is quite possible that 
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the discrepancy will simply disappear. If theory and experiment were to remain 
inconsistent, the discrepancy could in principle be an indication of new physics. 15 ) 

2.1.2. D — ► ttIi/, D — > Klu 

The shape of D semileptonic decay was predicted by lattice calculations 16 )' 17 ) 
before its precise measurement, and subsequently confirmed by the Focus, BaBar, 
and Belle experiments. Fig. [3] shows the shape of the form factor f+(q 2 ) for the 
decay D — > Klv calculated on the lattice 16 -*' 17 ) (orange and yellow bands), and as 
measured in experiment (points) by Belle, 18 ) BaBar, 19 ) and CLEO-c. 20 )' 21 ) 
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Fig. 3. The shape of the form factor f+(q 2 ) predicted with lattice QCD (orange and yellow bands), 
and as measured by Belle (green diamonds), by BaBar (magenta squares), and by CLEO (purple 
triangles) . 



The magnitude of the form factors can be used to determine \V c d\ and \ V CS \- At 
the conference Flavor Physics and CP Violation 2009, using CLEO-c data and the 
2005 lattice results, 16 ) Xin obtained 22 ) 

\V CS \ = 0.985(11)(103), (2-5) 
\V cd \ = 0.234(7) (25), (2-6) 

where the errors are from experiment and theory. Improved lattice calculations are 
underway for the theoretical errors. Prospects for the theoretical uncertainty to be 
significantly reduced are very good. 
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2.2. V cb and V ub 

These can be directly determined from the leptonic and exclusive semileptonic 
decays of B mesons. 

2.2.1. B — > Dlv , B — > 

In analyses that assume the unitarity of the CKM matrix, \V cb \ enters with a 
high power. For example, the kaon mixing parameter e is proportional to |Vd,| 4 , 
and \Vcb\ contributes as much uncertainty to the theoretical prediction for e as the 
kaon mixing parameter Bk, in spite of the fact that \V c b\ is known to much higher 
precision. 

Fortunately, it is possible to determine the form factors in B — > Dlv and B — > 
D*lv decays much more precisely than is possible in most leptonic and semileptonic 
decays. The case of B — > Dlv is especially simple to analyze. It is possible to connect 
\Vcb\ with the semileptonic form factors via a quantity in which the uncertainties 
cancel almost completely in the heavy-quark symmetry limit. The form factor may 
be obtained from the double ratio: 23 ) 



2 = < ^|c 70 6|i3 > < ffl 70 c|^ > 
< D\c 7o c\D > < /,'/,-„/, B > 



The factors in the denominator are used in the vector current renormalization re- 
quired in this amplitude. This ratio approaches one in the heavy-quark symmetry 
limit. In lattice calculations, uncertainties cancel almost completely in this limit. For 
physical values of the b and c quarks, uncertainties are proportional to deviations 
from the symmetry limit to high precision. h+ is the dominant term in the function 
Qb^d for the decay B — > Dlv. An unquenched lattice calculation gives 17 ) 

0b-i?(1) = 1.074(18)(16). (2-8) 

Using 1^*10(1) = (42.4 ± 1.6) x 10~ 3 from HFAG, 24 ) this produces 

\V A \ = (39.5(1.5) cxp (0.9) thco ) x 1(T 3 , (2-9) 

where the theory errors have been added in quadrature. The total error for |V^,| 
from B — ► Dlv is dominated by experiment. 

Experimental errors for B — > D*lv are smaller. From the Review of Particle 
Physics, 25 ) IKaIJ^I) = (35.9 ± 0.8) x 1CT 3 . Hashimoto et al. 26 ) defined a somewhat 
cumbersome combination of ratios from which .F(l) can be determined. Bernard et 
al. 27 ) investigated the much simpler ratio 

2 _ < D*\ci 3 ^b\B > < B\h jl5 c\D* > 
" <D*\c^c\D*> <B\h,b\B> (240) 

from which ^"(1) can also be determined. While calculational uncertainties do not 
cancel in this ratio as completely as in the ratios introduced earlier, they cancel to a 
high accuracy, and because this quantity can be calculated ten to twenty times faster 
than the set of ratios introduced earlier, 26 ) determinations of JF(1) from this ratio 
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are ultimately more accurate. Bernard et al. obtain = 0.921(13) sta t(20) sys . 

Combining with the experimental result yields 

\V cb \ = (38.7(0.9) cxp (1.0) theo ) x 10- 3 . (2-11) 

2.2.2. B -> ttIu 

Comparison between theory and experiment for B — ► iriv has been more trouble- 
some than for other lattice calculations in CKM physics. Leptonic decays and BB 
mixing amplitudes are described by a single parameter. The semileptonic decays 
B -> D^Hu and K -> ixlv can be described to high accuracy by a normalization 
and a slope. For B — > n£v, on the other hand, the form factors have a complicated q 2 
dependence. Lattice data have covered only the low momentum, high q 2 end of the 
pion momentum spectrum, and errors are highly q 2 dependent and highly correlated 
in both theory and experiment. 

It has long been understood that analyticity and unitarity can be used to con- 
strain the possible shapes of form factors. Consider mapping the variable q 2 onto a 
new variable, z, in the following way: 

where t + = (rrifl + m^) 2 , t_ = (mg — m n ) 2 , and to is a f ree parameter. Although 
this mapping appears complicated, it actually has a simple interpretation in terms 
of q 2 ; this transformation maps q 2 > t + (the production region) onto \z\ = 1 and 
maps q 2 < t+ (which includes the semileptonic region) onto real z G [—1,1]. In the 
case of B — > nils, the physical decay region is mapped into roughly —0.3 < z < 0.3. 
In terms of z, the form factors can be written in a simple form: 

1 00 

Most of the q 2 dependence is contained in the first two, perturbatively calculable, 
factors. The Blaschke factor P(q 2 ) is a function that contains subthreshold poles 
and the outer function cj)(q 2 ,to) is an arbitrary analytic function (outside the cut 
from t+ < q 2 < 00) which is chosen to give the series coefficients a simple form. 
See 28 ) ; 29 ) and references therein for the explicit forms of these expressions. With 
the proper choice of (p(q 2 ,to), analyticity and unitarity require the a/% to satisfy 



N 



ai<~\. (2-14) 

k=0 



The fact that —0.3 < z < 0.3 means that according to analyticity and unitarity, only 
a few terms are required to describe the form factors to 1% accuracy. 

Calculations have been performed by the Fermilab Lattice and MILC collabora- 
tions 28 ^ using Fermilab b quarks, and by the HPQCD collaboration using NRQCD 
b quarks. 30 ) Figure [4] shows the result from Fermilab/MILC of a fully correlated 
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Fig. 4. Results for the normalized B — > nlv form factor P+(f>+f+ from the Fermilab/MILC lattice 
calculations (circles) and BaBar (stars). The solid red line is the results of a fully correlated 
simultaneous fit. Requiring that lattice and experiment have the same normalization yields 

\V ub \. 



simultaneous z-&t to the Fermilab/MILC lattice data and the BaBar 12-bin experi- 
mental results, 31 ^ with \V u b\ being a parameter in the fit. Good fits may be obtained 
with just three terms in the expansion to the lattice or the BaBar data separately, 
and to the combined lattice and experimental data, with the result 

\V ub \ x 10 3 = 3.38 ±0.36. (2-15) 

Because the 10% uncertainty comes from a simultaneous fit of the lattice and ex- 
perimental data, it contains both the experimental and theoretical errors in a way 
that is not simple to disentangle. If we make the assumption that the error in \V u b\ 
is dominated by the most precisely determined lattice point, we can estimate that 
the contributions are roughly equally divided as ~ 6% lattice statistical and chiral 
extrapolation (combined), ~ 6% lattice systematic, and ~ 6% experimental. The 
largest lattice systematic uncertainties are heavy quark discretization, the perturba- 
tive correction, and the uncertainty in g B * Bn , all about 3%. 

HPQCD has obtained compatible results using NRQCD b quarks and asqtad 
light quarks. 30 ^ They obtain consistent results from several types of fitting methods. 
Applying their results to 2008 data from HFAG 32 ) yields 

\V ub \ x 10 3 = 3.40 ± 0.2012:11, (2-16) 

where the first error is experimental and the second is from the lattice calculations. 
2.2.3. f B &ndf Bs 

f B , and f Bs have been calculated by HPQCD 33 )- 35 ) and by Fermilab/MILC. 10 )' 11 ) 
HPQCD obtained f B = 190(13) MeV and f Bs = 231(15) MeV, and Fermilab/MILC 
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obtained f B = 195(11) MeV and f Bs = 243(11) MeV. It is complicated to know 
how to combine the uncertainties since the two calculations share some uncertain- 
ties. One could simply take the smaller of the two individual uncertainties as the 
combined uncertainty and use for an average f B = 193(11) MeV and f Bg = 237(11) 
MeV. 

The direct observation of the leptonic decay of the B into a r and its neutrino 
has become increasingly precise. The world average 25 ^ of r(j + u) //total = 1-4(4) x 
10 -4 allows a determination of \V U \,\ to 14%. While not (yet) competitive with the 
determination for B — > irlv, this is a level that would have been hard to imagine a 
few years ago. 

2.3. V ud and V us 

\V u d\ and \V US \ are directly related to the physics of pions and kaons. The top 
row of the CKM matrix provides a precise test of the unitarity of the matrix via the 
relation 

\V ud \ 2 + \V US \ 2 + \V ub \ 2 = 1. (2-17) 

Since \V u b\ is around 4 x 10~ 3 , the third term's effects are negligible in this relation. 
\V U( i\ can be obtained from pion leptonic decay and lattice calculations of the pion 
decay constant, but these determinations are not as accurate as those from nuclear 
beta decays. Lattice calculations of f n provide tests of lattice techniques. From 
nuclear beta decay, we have \V u d\ = 0. 97425(22). 36 ) This contributes an uncertainty 



of ±0.00043 to Eq. (2-17). Hence, a comparable level of uncertainty is of interest in 
\V US \, an uncertainty of around 0.5%. Remarkably, because \ V US \ may be determined 
from ratios in which many of the theoretical uncertainties cancel and the remain- 
ing uncertainties are proportional to SU(3) flavor breaking, this daunting level of 
accuracy is not out of reach. 

\V US \ may be obtained in two ways from ratios. |Vu S /VW| may be obtained via 
the ratio of pion and kaon decay constants fn/fx- \Vus\ may be directly obtained 
from a double ratio similar to that used to obtain \ V c b\. 

2.3.1. /,,■ and f K 

In 2004, Marciano emphasized the possibility of obtaining (V^/V^l from the 
ratio of the leptonic decay constants. 37 ^ He provided a necessary radiative correction 
to the ratio and obtained 

^4^ = 0.2757(7). (2-18) 

\Vud\ Jtt 

At Lattice 2008, kaon physics was reviewed by Lellouch. 5 ) Most lattice uncer- 
tainties cancel in the ratio fx/ fir- The dominant uncertainty results from a combi- 
nation of chiral extrapolation and statistical error. Lellouch performed a weighted 
average of the unquenched results and obtained fx/ fw = 1.194(10). When combined 
with Marciano 's result, this yields 

0.2309(18), (2-19) 



\V ud \ 

or about 0.8% uncertainty, dominated by theory. 
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2.3.2. K^itlv 

Like B — > Z?/z^, the semileptonic form factor for the decay -ftT — > 7r/^ may be 
obtained from a double ratio. In the case of K — ► 7rZi/, flavor SU(3) symmetry 
ensures that most statistical and systematic errors cancel in the symmetry limit, so 
that errors are very small in the physical case. One starts from the double ratio 



, 2^MX <*\V \K><K\V Q \*> 

JOWmaxJ (M K + M v )<*\V \*><K\Vo\K>' 1 ' 

The denominator factors remove current renormalizations. As shown by Becerivic 
et al., 38 ) by interpolating in q 2 to obtain fo{q 2 = 0) and using /+(0) = /o(0), one 
obtains the required form factor /+ to high accuracy. Lattice results for /+(0) were 
reviewed in Lellouch, 5 ) with the average /+(0) = 0.964(5). The experimental result 
for K13 decay is 39 ) \V US \ x /+(0) = 0.21664(4). This gives 

\V US \ = 0.2247(12), (2-21) 

or about 0.5% accuracy. 

Thus, the accuracies obtained for the unitarity test of Eqn. (2-17) via these two 
methods are reaching the required 0.5% level. Because the quantities involved are 
unusually simple lattice quantities, and because the necessary information can be 
obtained from ratios in which most uncertainties cancel, the prospects are good for 
surpassing the required accuracy soon. 

In 1984, Leutwyler and Roos made the estimate 40 ^ /+(0) = 0.961(8). This es- 
timate has stood the test of time remarkably well. It is in agreement with modern 
lattice calculations, and lattice calculations have only recently surpassed it in preci- 
sion. Some lattice calculations, such as heavy meson decay constants, have revealed 
that QCD predicts something quite different from the old conventional wisdom that 
had been arrived at via quark models. In K — > ttIu decay on the other hand, 
Leutwyler and Roos's rough estimate of twenty years ago has proved remarkably 
robust in the face of more fundamental calculations. 

2.4. V t d and V ts 

Vtd and Vts cannot be determined from tree level processes, but can be deter- 
mined from loop effects in meson-antimeson mixing. BB and B S B S mixing depend 
on Vtd and Vts in a simple way. The kaon mixing parameter e depends on a compli- 
cated mixture of quantities including Vtd and Vt s , and I will discuss it in this section, 
too. 

2.4.1. BB and B S B S mixing 



In jB( a) mixing, the quantities f B \[B~ B , f Bs yj B Bs , and £ = f Bs y/ B B J (/bVBb) 
are calculated. The mixing amplitudes likef B \/B B are traditionally parametrized by 
a decay constant and a "bag" parameter, but the amplitude is actually calculated as 
a single quantity. Many of the usual lattice uncertainties cancel to a large extent in £. 
Its uncertainty is dominated by a combination of chiral extrapolation and statistics. 
f Bs y B Bs has only a mild chiral extrapolation (arising from the sea quarks), but it 
has a larger share of the usual lattice uncertainties, such as discretization, operator 
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matching, and lattice scale. £ and fs s \/Bb s therefore have relatively uncorrelated 
uncertainty budgets, and are the most suitable pair of quantities for inputs into 
global fits. Jb^/Bb, by contrast, has a larger admixture of all of these uncertainties. 
The expressions for B^ mixing can be written 

AM d(s) oc \V: d{s) V tb \ 2 fl d(s) B Bd{aV (2-22) 
where < B° s \O s L {d) \B° s >= IM^I^Bb^), and 

O s £ d) = [5*7^(1 - 7 5 )s i (^)][6 i 7^ 1 (l - Ts)^(dJ)]. HPQCD have published results for 
the mixing amplitudes 33 ) f Bs ^/B~B~ a = 266(18) MeV and f B V^B = 216(15) MeV. 

The most precise constraint on the p-rj plane from B and B s mixing can be 
obtained from the combination 



\V U \ y AM S Ms d V ' 

For £, HPQCD has 33 ) £ = 1.258(33). A preliminary result from Fermilab/MILC gives 
£ = 1. 205(50). 41 ) The uncertainty in this second number is expected to drop in the 



final result expected later this year. Combining the HPQCD result with Eq. (2-23) 
gives 

||4 = 0.214(1)(5), (2-24) 
I Vts\ 

where the first error is experimental and the second is theoretical. In spite of the high 
precision already achieved theoretically for £, the even higher experimental precision 
for AM d j AM S means that a much larger payoff in constraining the CKM matrix 
will result from further improved lattice calculations, which should be possible. 

One of the most important possible outcomes of high precision CKM constraints 
would be an inconsistency in the constraints, indicating new physics. The fact that 
B and B s mixing originate in loop effects has led to speculation that they may be 
more sensitive to new physics in the loops than are tree-level processes. In BB 
mixing, there are five operators in the most general basis for four-quark scattering 
that can arise in BSM physics. These have been tabulated and evaluated in the 
quenched approximation by Becirevic et al. 42 ) It is straightforward to repeat these 
calculations in unquenched calculations, and efforts to do so are underway. 

2.4.2. A comment on DD mixing on the lattice 

The calculation of the short-distance part of DD mixing is identical to the BB 
mixing case just discussed. These calculations have received less attention on the 
lattice than the analogous calculations in the B system because in DD mixing, there 
are long-distance contributions from intermediate states. To calculate these on the 
lattice in a naive way would require a four-point function, as opposed to the three- 
point functions used in current BB mixing calculations. This would require a number 
of numerical operations proportional to the lattice volume, V, squared, as opposed 
to the 0{V) operations used in current calculations. This is beyond the capability of 
current computers. It is probably possible to attack long-distance effects with more 
intelligent methods, but these have not been worked out. 
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2.4.3. KK mixing 

The formula for the KK mixing parameter e is 




+ 7 ? ^(x t )i m [(^ s ^ 2 ] + 2r ?ct s(x C) x t )i m [(y cs y; d y ts v t * d )]]. 



(2-25) 



This depends on the renormalization group invariant bag parameter Bk and on the 
CKM matrix elements V cs , V c d, Vt s , and Vtd and on several other correction factors 
with uncertainties. Many of these uncertainties contribute to the uncertainty in the 
bound that kaon mixing can give in the p-r/ plane, which makes the prospects for 
improving the bound more problematic to estimate. 

World results for Bk were reviewed by Lellouch at Lattice 2009, who obtained 5 - 1 ' 43 ) 
Bk = 0.723(37). The largest uncertainties in constraining the p-fj from e arise from 
B K and V cb . 25 ^ (The latter is in the sense that, assuming CKM unitarity, (Vt s V^) 2 
is proportional to the Wolfenstein parameter A 4 , which in turn is proportional to 
|Vrf,| 4 .) However, as illustrated in Lellouch, 5 ) even if these uncertainties were to be 
reduced to zero in expression for e, there would still be significant uncertainty from 
the rest of the parameters. In Fig. [TJ one sees that the constraint in the p-rj plane 
from e is less stringent that those from B and B s mixing and from V u b inclusive and 
exclusive. Improving this bound will be possible, but complicated. 



This is the only element of the CKM Matrix about which lattice QCD has 
nothing to say. 



3.1. The current generation of lattice QCD calculations 

The current generation of lattice QCD calculations serves both to demonstrate 
the correctness of lattice calculations and to deliver quantities for physics that can 
only be obtained with lattice QCD. Foundations are also currently being laid for cal- 
culations in strongly interacting Beyond-the-Standard-Model theories, should they 
be discovered. One obvious goal of the current generation of lattice calculations is to 
demonstrate the correctness of lattice methods, such as by correctly post-dieting the 
hadron spectrum. Some of the weak matrix elements I have discussed can also serve 
to do this by constructing CKM-in dependent ratios of physical quantities. Occasion- 
ally, it has been possible to make predictions, even in low energy hadronic physics, 
such as in the shape of the form factor in D s — * Klv decay that I discussed. Also 
of greater interest to the general community may be the resolutions of the contra- 
dictions that have occurred in a few places between lattice results and other results. 
The lattice results for fjj s lie significantly below the experimental results. There are 
also interesting two-sigma tensions between inclusive and exclusive results for \V u b\ 
and \ Vcb\- It is not known at present how these will play out. 

The lattice determinations of CKM matrix elements that are the topic of this 
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article are among the most important current deliverables of lattice gauge theory for 
particle physics and for the world of physics in general. Some of the the fundamental 
parameters of the Standard Model can be directly determined only with lattice QCD, 
including the light quark masses and the constraints on the CKM matrix from KK, 
BB, and B S B S mixing. 

Some flavor issues may have their resolution in Beyond-the-Standard-Model 
physics. Lattice calculations are beginning to address these. One example that 
has become tractable with the current technology of lattice QCD is the search for 
near-conformal behavior in gauge theories. 44 ) One explanation for the smallness of 
flavor-changing neutral currents in the context of Technicolor models of dynamical 
breaking of the weak gauge symmetry involves an almost conformal fixed point in 
the symmetry breaking sector. Lattice calculations have begun to search for such 
theories. 

3.2. The coming generation of lattice calculations. 

I have not said much about prospects for future improvements in the results I 
have discussed. This is because the prospects for almost all of them are excellent. 
Most of the most important lattice results in CKM physics involve hadronically sta- 
ble mesons. These are the simplest hadrons to study on the lattice (or elsewhere), 
and there are no known impediments to higher precision. Improvements to meth- 
ods, algorithms, and computers are putting errors such as discretization and chiral 
extrapolation under increasingly good control. Uncertainties due to operator nor- 
malizations fall more slowly with brute force computer power, but here also we are 
in good shape. Lattice perturbation theory converges more or less as well does di- 
mensionally regulated perturbative QCD. 45 ) On the lattice, however, there are also 
nonperturbative renormalization methods available that are not available in the con- 
tinuum, 46 ) putting short-distance QCD under better control on the lattice than in 
the continuum. Achieving the sub-per cent precision required by experiment in KK, 
BB, and B S B S mixing for bounding the ~p-fj plane is by no means an unreachable 
goal. 

An increasing variety of physical quantities will become tractable in the coming 
years. I have focused in this article on processes with single mesons in the final state, 
because these are the most precise lattice CKM calculations. Exclusive processes 
with more than one hadron are more complicated because the translation between 
Euclidean final states and Minkowskian final states are more complicated than for 
single-hadron states, but methods exist for some of these. First-principles methods 
have been worked out for the two-pion decays of K mesons, which are needed to 
extract CKM properties from e'/e. 47 -* Chiral perturbation theory can be used to 
help analyze final states, since the decay pions have relatively small momentum. 
These calculations are somewhat more cumbersome than single hadron calculations 
and will require more computation, but can ultimately be done with current meth- 
ods, and prospects for progress are good. For interesting decays like B — > irn, on 
the other hand, no practical methods exist, and new methods must be developed. 
The long-distance contributions to DD mixing have not yet been addressed with 
lattice calculations. The well-advanced topics I discussed in the previous section are 
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done with two- and three-point correlation functions. The higher-point correlation 
functions required in DD mixing and similarly in the "light-by-light" scattering non- 
perturbative contributions to g — 2 for the muon require more advanced methods. 
Those for light-by-light scattering are now being worked out. 

Some Beyond-the-Standard-Model strongly coupled gauge theories are within 
the reach of current methods, as already discussed. 44 ) Others, such as the simplest 
strongly coupled supersymmetric gauge theories are tractable with straightforward 
extensions of current methods and are being actively investigated. 48 - 1 Others will re- 
quire more significant improvements to methods. These are all being given increased 
attention as the LHC era dawns. How new data from the LHC will affect these de- 
velopments is impossible to foresee. With increasingly precise results and a broader 
variety of applications, lattice gauge theory will continue to play an expanding role 
in particle physics in the LHC era. 
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